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Abstract 

This paper explores finding the number Uh of undirected hamiltonian 
paths in an undirected graph G = {V,E) using lumped/ideal circuits, 
specifically low-pass filters. Ideal analog computation allows one to com¬ 
puter Uh in a short period of time, but in practice, precision problems 
disturb this ideal nature. A digital/algorithmic approach is proposed, and 
then it is shown that the approach/method operates under theoretically 
feasible (polynomial) time: 0(n^®“) + 0{ndPC{p2,nd)), where n = \V\ 
and PC{p 2 , Ud) refers to the time complexity of finding a root of a polyno¬ 
mial with polynomial coefficients being p 2 digits (with coefficients limited 
within 2^^^^ and in magnitude) and Ud referring to number of 

degree, with p 2 = and Ud = 


1 Introduction 

This paper is all about converting an undirected graph G = {V, E) to a signal 
processing problem, whether that be using a digital or analog computer, and 
use the reduction to resolve a ^P-complete problem - here the counting version 
of the hamiltonian path problem - efficiently/feasibly, at least if we interpret 
“feasible” as “polynomial time relative to \V\ = n.” 

Graphs and functions have been studied together - for example, quantum graphs. 
In this line, this paper is one of these cases. However, it is surprisingly difficult to 
find an example where a graph is represented as a function and signal processing 
tools are used to analyze the function. 

Many methods used in this paper are standard analog signal processing tools, 
and if there is anything novel in this paper, it would be about adapting the 
analog signal processing method to meet the numerical needs of the problem 
this paper intends to tackle. 

Thus, the main reference for this paper is essentially analog signal processing 
books, or possibly even textbooks. I assume that the readers are familiar with 
signal processing terminology - terms like HR, FIR, steady-state and transient 
response and so on. Or more precisely, the readers are assumed to be familiar 
with sinusoidal signal processing analysis. 
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2 Algorithmic approach 

In this paper, all graphs are assumed to be undirected. Key primary insights are 
written in bold typeface, and secondary insights that are nevertheless of some 
importance will often be listed in itemized lists. 

2.1 Two Goals 

The first goal of this paper is to show that the counting variant of the hamilto- 
nian path problem is effectively a signal processing problem. By reduction, this 
means that several important computation problems are “isomorphic” to signal 
processing problems. 

The second goal of this paper is exploring whether feasible computation of the 
number of hamiltonian paths in an undirected graph is possible, via signal pro¬ 
cessing understanding of the hamiltonian path problem. For analog computa¬ 
tion, the concept of “feasible” is obscure, but for digital computation, “feasible” 
will refer to computation in polynomial time. 

For both cases, the answer is yes. 

2.2 Two Core Concepts 

These two core concepts are used to convert a hamiltonian path counting prob¬ 
lem into a signal processing problem. The summary of core concepts is provided 
at the end of the next subsection, as part of Conclusion. 

2.2.1 Core Concept 1 (CCl): reducing a graph to a function f(t) 

As the title of this sub-subsection says, the first core concept is reducing an 
undirected graph to /(t) : M —>■ C. 

Initially, one forms x(t) before forming f(t). Each vertex is assigned a sep¬ 
arate angular frequency. Then each walk of a graph restricted to containing 
\V\= n vertex visits is also assigned an angular frequency. By the definition of 
a walk, visiting one vertex more than once is allowed. Let us call such walks as 
n-walks for simplihcation. 

Each n-walk can be assigned angular frequency by simply summing up angular 
frequency of each vertex by visit order. Thus, the angular frequency of a n-walk, 
expanded as the sum of vertex terms, may have one vertex added more than 
once. 

The main reason for this function conversion is that in this process, 
all hamiltonian paths share the same angular frequency, because they 
have to contain all vertices in an undirected graph. Thus, one would like 
to have walks that are not hamiltonian paths to not share the angular frequency 
of hamiltonian paths. This is done by using particular angular frequency as¬ 
signing methods so that each n-walk that shares the same set of vertices along 
with the same vertex-corresponding visit frequencies, has the same angular fre¬ 
quency. Of course this is not the only way, but this is the most straightforward 
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way of proceeding. 

By visit frequency, if a n-walk visited some vertex Va twice, then the visit fre¬ 
quency of Vo is 2. 

Then I “shift” angular frequencies so that hamiltonian path angular fre¬ 
quency is zero, which forms y(t). This allows one to formulate a hamil¬ 
tonian path counting problem as the signal processing problem of 
finding the period-adjusted average - or simply, the zero frequency 
amplitude - of the Fourier series y(t). 

For convenience, I scale time/angular frequency, which gives one a Fourier series 
form of an input graph: /(<). 

2.2.2 CCl: more detailed conceptual explanations (CCIDI) 

Definition 2.1 (visit frequency/occurrence). Vertex v visit frequency/occurrence 
of a walk refers to the number of times v is visited in a walk. Visiting order 
does not matter. 

Definition 2.2 (n-walk). A n-walk is a walk with n vertex visits. That is for 
G = (V, E), if one sums up visit occurrence of all vertices of a walk, one gets n. 

Definition 2.3 (visit pair). Visit pair for vertex v in a walk is the pair {v,k) 
that shows visit occurrence by natural number k. Visit pair of a walk is defined 
as the set of all such visit pairs, each corresponding to each vertex, for graph 
G. For example, for a graph with 4 vertices with V = {1,2,3,4}, a walk may 
have visit pair of {(1, 2), (2,0), (3,1), (4,1)}, where 4 = 2-|-0-|-l-l-l. 

A walk is assumed to be a valid walk allowed by undirected edge construction 
of G. The definition for a walk follows standard graph theory terminology. 

The first motivation for constructing a grid representation f{t) for a graph comes 
from the following list: 

• A set of n-walks that share the same visit pair is assigned same angular 
frequency. That is, each walk is assigned angular frequency, but walks 
with a common visit pair have same angular frequency. 

• Each n-walk has amplitude of 1. That is, each n-walk with angular fre¬ 
quency oj is said to contribute e“h 

• One wishes to ensure that a different n-walk visit pair is assigned a sepa¬ 
rate frequency, distinguishable from others. 

• If the second item in this list is ensured, then it is possible to distinguish 
a n-walk visit pair that captures hamiltonian paths, where all n vertices 
of graph G exist, each with visit occurrence of 1. 

• For hnal graph grid representation f{t), we would like to make zero fre¬ 
quency hamiltonian path frequency. This is to extract the number of 
hamiltonian paths rih using low-pass filtering techniques. 
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For this section, let U be the set of angular frequencies with corresponding n- 
walk in a graph. 

Now consider how one assigns angular frequency to each walk. 

Definition 2.4 (Vertex-number). A vertex-number is assigned to each vertex. 

If not mentioned explicitly, vertex refers to a vertex-number it is assigned. In 
a particular implementation of the method used in this paper, each vertex is 
V = n* with i G 'L+. Thus, V = {n, n^, n ^,.., n”}. 

From now on, V = {n,n?,n^, ..,n”} is used. 

Definition 2.5 (Assigning angular frequency to a n-walk). Following from the 
visit pair p of a walk, with V in mind, recall the definition of p: p = {(u, k)\v G 
V}, where k is visit occurrence. The angular frequency oj of each path is a; = 
kv. To state again, each walk contributes e“b 

Definition 2.6 (A„). For this sub-subsection only, A„ refers to the car¬ 
dinality of the set of all the n-walks sharing angular frequency u. Separate 
definitions will be provided whenever is used in different contexts. 

Definition 2.7 (x(t)). x(t) be the sum of all n-walks. 

That is, x(t) = 

Definition 2.8 (a^,). That is, refers to hamiltonian angular 

frequency for x{t). For y{t) and f{t) this will not be the case. 

Definition 2.9 {y{t)). yit) = a;(t)e“*“'*‘. This shifts all angular frequencies of 
x{t) by ah, and hamiltonian angular frequency for y{t) is now zero. 

Definition 2.10 (maximum angular frequency). For this definition only, I 

will relax the definition of A„ and U. U is assumed to be hnite, and contains 
real numbers. Let A„ be some complex number that depends on u. Let some 
signal k{t) = Then maximum angular frequency of k(t) refers to 

|u| (notice the absolute sign, and this is not mistake) with greatest 1^1 in the set 
U. The definition for maximum angular frequency applies for all parts of this 
paper. 

Definition 2.11 (signal/part of angular frequency ui). Consider /x(t) = + 

f^, ‘ duj'. Then signal/part of angular frequency uj of y(t) is {A^ + duj) 

Definition 2.12 (amplitude at w of some function g{t)). It is A^ + B^div 
mentioned in the definition for signal/part of angular frequency oj (right above). 

Definition 2.13 (,/(t)). f{t) = y{ct). c in this implementation is defined as 
c = \ Thus /(t) will have maximum angular frequency of n^" )+n+i^ 

with minimum non-zero angular frequency (with non-zero amplitude) being 
n^" K [In practice, maximum angular frequency of f{t) will be lower, and 
this is over-estimate. Similarly, minimum non-zero angular frequency of f(t) is 
under-estimate.] 



Signal Processing Understanding of Hamiltonian Paths 


5 


This completes defining process of f(t), which is to be used for low-pass 
filtering to obtain Uh- However, actual implementation of how to obtain f(t) 
may be inferred but has not been presented yet. That is, we know that f(t) is 
essentially the sum of different walks, but how will different walks be combined 
efficiently, so that we do not have to check individual walk to get f(t) for some 
tl 

2.2.3 CCl continuation: efficient walk sum 

For now, I will give analogy to a digital circuit that can easily be understood 
back in terms of algorithms. The idea is first given for x{t), as calculation of 
y{t) and f{t) essentially follow automatically from calculation of x{t). 

• Recall that vertices were assigned vertex-numbers. 

• Each vertex has two storage parts: incoming and outgoing storage. 

• The idea now is to treat the entire graph as a digital synchronized/timed 
circuit. 

• The digital circuit is assumed to be edge-triggered. 

• Wires between vertices are connected to edge connectivity. If (^ 1 ,^ 2 ) G E, 
then a direct wire connecting vi to V 2 exists. Note that we assume two 
wires exist for {vi,V 2 )' one wire wi sends signal from vi to V 2 , while the 
other wire W 2 sends signal from V 2 to vi. 

• wi is incoming wire to V 2 , outgoing wire to vi. W 2 is incoming wire to vi , 
outgoing wire to V 2 - 

• When the digital circuit is triggered, each vertex first copies the signal in 
its incoming storage to outgoing storage. Then each vertex activates all 
the wires previously deactivated, so that signal may be transmitted. 

• When the wires are activated, each vertex v receives signals from incoming 
wires, sum up all signals, multiply the sum by and store the result to 
incoming storage. 

• Right after the circuit is edge-triggered, each vertex also transmits signal 
stored in outgoing storage to outgoing wires. 

• After new data are copied to incoming storage, all wires and the circuit 
are deactivated, waiting or the next trigger to occur. 

• Before the first trigger, initial data stored in the incoming storage of each 
vertex v is e*"*. 

• After n — 1th trigger time passes, one simply sums up incoming storages 
of all vertices and outputs x{t). 
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• As the last point, note that data in the storages will be vectors of poly¬ 
nomial coefficients that are obtained via Taylor-expanding around t = 0. 
Refer to the sub-subsection “The polynomial form of f(t)” for more de¬ 
tails. 

(At this point, one may be reminded of Boltzmann machines, though I will not 
go over this.) 

Essentially, the idea above uses the following characteristic: 

(g™l,l‘ _|_ g™l,2t^g™2,3t _ g*(*'l,l+'“2,3)* _|_ gi(l’l,2-|-l'2,3)t 

I'l l and Vi ^2 can be thought of the first vertex visited for each corresponding 
n-walk. Both walks share the same second vertex visited, which is V 2 , 3 - Notice 
that the first index (1 in vi_ 2 , for example) presents the visit order of a walk, 
while the second index distinguishes actual vertex. Thus, and ^ 3^3 represent 
the same vertex, but with a different visit order. 

The full details will be given, but for now let us illustrate the principles using 
the example in Figure and 


Figure 1: A 4-vertex graph 



In Figurej^ because the original graph in Figure[l]is a 4-vertex graph, there 
are four layers or four depths, labelled with L1,L2,L3,L4. Each layer i contains 
all vertices in a graph. A vertex v in layer i is connected to a vertex w in layer 
i1 whenever {v,w) € E. This grid procedure effectively simulates an actual 
walk. 

At Layer 1 (LI), each vertex v transmits e™* to the edges, or wires, it is con¬ 
nected with. These edges connect to the vertices at the next layer L2. For other 
layers, each vertex v sums up all the function it received from the wires starting 
from the previous layer and multiplies the sum by e*"*, and then transmits the 
product to the wires that connect v with the vertices of the next layer. The 
final layer L4 has an additional step, since there are no wires that connect to 
the next layer in Figure]^ Instead, all results obtained at each vertex at L4 is 
summed up, which results in x(t). 

To summarize, the vertices in LI always act as oscillators, every edge that con¬ 
nects one layer to the next layer acts as a right-directional wire without any 
transmission delay and the vertices in each layer except LI work first as a sum¬ 
mer and then a multiplier coupled with an oscillator. After the final layer, a 
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LI 


L2 


L3 


L4 


summer adds up all results in the vertices in the final layer. 

The constructed x{t) does not have zero frequency as a hamiltonian path fre¬ 
quency and thus multiplication by where ah is the hamiltonian path 

angular frequency of x{t), is needed to produce y(t). Then for convenience an¬ 
gular frequencies may be scaled by multiplicative factor (just by changing time 
scale) to produce f{t). 

This completes introduction of the grid, and in fact, this is all there is for the 
grid. 

2.2.4 CCl implementation: full details 

Definition 2.14 (2'+, Z~, Z-). or Z+ refers to the set of positive 
integers. Similarly, Z~ refers to the set of negative integers. 

Definition 2.15 (“less than”, “more than”, “greater than”, “smaller than”). 
Unless otherwise noted, these are all comparisons in magnitude/size/absolute 
value. 

Definition 2.16 (Base-n expansion). Base-n expansion of some number k is 
basically expressing k in base-n: fc = ± with 0 < Op < n. 

Definition 2.17 (graph, n). A graph G is denoted with G = (U, E) as done in 
the standard literature, n = |U| is assumed whenever n appears. 
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Definition 2.18 (walk, n-walk, hamiltonian path). A walk is defined as in the 
standard graph theory vocabulary. A walk that has n vertices is called n-walk. 
Let us represent a walk with a list (tuple) of vertices in a traversing order from 
the start vertex to the end vertex. By the definition of a walk, one vertex can 
appear more than once in a list. A hamiltonian path, as defined in the standard 
graph theory vocabulary, is a walk with n distinct vertices, where \V\ = n. 

Definition 2.19 (vertex). A vertex is assigned a number. Each distinct vertex 
has a distinct vertex-number. Let V = {n, n^, n^,..., n"}. From now on, one 
can assume a vertex as a number whenever appropriate. 

Definition 2.20 (n/i, Up). nt is the number of hamiltonian paths of G. Up is 
the total number of n-walks of G. 

Definition 2.21 (Walk-number). The walk-number of a walk is defined as the 
sum of all elements (vertices) in the list of a walk. Note that “list” is used to 
refer to the fact a same vertex may be visited several times and it may need to 
be summed up several times. 

Recall that the walk-number of a walk represents the angular frequency of 
a walk in x{t). It is certainly possible that two walks may occupy the same 
frequency. If there are k walks that occupy the same frequency uja, then the 
amplitude at the frequency would be k in Fourier series language, or k6{u! — uJa) 
in Fourier transform language where 6{(jj) is a dirac delta function. 

The maximum number of vertices inside a walk is restricted to n, for 
sake of convenience. 

Lemma 2.1. Given V as defined above, each walk-number has a unique visit 
pair attached. That is, no other visit pairs may generate that walk-number. 
(That each visit pair has a unique walk-number is trivial to prove.) 

Proof. The proof is simply the basis representation theorem, where basis are 
elements in V. One exception to this proof, though, arises when a list f repre¬ 
senting a walk may be of {k, k, ...k) with |^| = n and k = n*, or in words, there 
are n k’s in f. In this case, nk = meaning the vertex-number f equals 

one of vertices in V. But this should not matter whenever walks one deals with 
have same number of vertices. □ 

Following from above: 

Definition 2.22 (Contribution of each n-walk to x{t)). From above, each walk 
has a walk-number k. Each n-walk is said to contribute e®** to x{t). 

Definition 2.23 (Amplitude). For any arbitrary function a{t) expressible as 
aft) = J2^=-oo where d is constant and does not vary with ui, is 

said to be amplitude of aft) at angular frequency uj. 
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2.2.5 CCl continuation: grid: x(t) 

Definition 2.24 (Grid, wires). A grid consists of n depths, with each depth 
being equivalent to a column. Each depth contains n vertices as in V. Each 
wire connects a vertex Va from fth depth to a vertex point of Vjs in z + 1th depth. 
A wire is connected between Va to vp if and only if {va,vp) G E. 

Definition 2.25 (Eunction transmission: first depth case). In the first depth 
(first column), each vertex Va transmits 

Definition 2.26 (Eunction transmission except for first and nth depth). Defin¬ 
ing for each Va in arbitrary zth depth. All incoming wire transmissions zc^(t) 
from each wire C from z — 1th depth to Va in zth depth are summed, or equiva¬ 
lently w\ = then multiply by and transmit Uy^ = to 

each wire starting from Va- 

Definition 2.27 (Vertex point function transmission: nth depth case). All 
incoming wire transmissions zc^(t) from each wire C from n — 1th depth to Va 
in nth depth are summed, or equivalently wx = wq. And then multiply by 
e““*, resulting in = e™°‘*'wx. Xideai{t) = is the output of the grid, 

not considering polynomial coefficient errors and polynomial degree truncation 
involved. 

Eor each depth z, shows the sum of all vertex-numbers representing 

z-walk. 

2.2.6 CCl continuation: post-grid: y{t) 

Simply, this post-grid procedure is all about calculating y{t) = a;(t)e“*“'** where 
Uh = X]r=i hamiltonian frequency of x{t). Thus, y{t) has 0 has hamilto- 

nian frequency. 

2.2.7 CCl continuation: post-grid: f{t) 
f{t) is defined as f{t) = y{ct). c was defined as c = 

Let the angular frequencies of f{t) be labelled with oj. uj = 0 refers to hamilto¬ 
nian frequency. 

2.2.8 The polynomial form of f{t) 

The polynomial form of y{t) is easily obtained by expanding e*"* of each vertex 
V to the imposed polynomial degree rzd.i (defined below) and do multiplication 
of the polynomial form of e*’'* (taylor expansion around t = 0) with the sum 
of incoming polynomials. The polynomial form obtained, truncated to degree 
Ud,It is passed to outgoing wires. This gives us final accuracy up to polynomial 
degree Ud^i- There is coefficient error involved during the calculations, and this 
will be dealt separately. 
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2.2.9 Core Concept 2 (CC2): Taking advantage of known informa¬ 
tion: signal processing tools, multi-step 

Definition 2.28 (polynomial term, polynomial form). The polynomial term 
of degree d of some function g(t) is the Adt‘^/dl term that exists in the Taylor 
expansion of g(t) around t = 0. The polynomial form of some function g(t) of 
the assigned degree Ud refers to the truncated Taylor expansion of g(t) around 
t = 0, where polynomial terms of degree greater than rid are thrown out. 

Definition 2.29 (term). The term of some function will refer to the polynomial 
term with other terms also assumed to be polynomial terms, if not explicitly 
stated otherwise. When one says term,” however, “term” may refer 

to transient response terms, and meaning will easily be inferred within contexts 
the word “term” is used. For example, the sentence like “A2t^e“*/2! term of 
function (/(t)” may assume, depending on context, the decomposition/dissection 
of g{t) into g{t) = la other contexts, g{t) may additionally 

include steady-state response, depending on contexts for the example sentence. 

Definition 2.30 (A). A does not really have a fixed meaning, and is used 
to represent an arbitrary number, amplitude and so on. Similarly with other 
letters without explicit definition like B, C, et cetera. 

Definition 2.31 (function, t). All functions are assumed to be of K —)■ C, with 
t being an independent variable. 

Definition 2.32 (finite-degree truncation error). Finite-degree truncation er¬ 
rors of the polynomial form of some g(t) are the errors caused by Taylor- 
expanding g{t) around t = 0 and truncating polynomial terms of degree greater 
than rid (or rid^ in case of step 1) (defined below). 

Definition 2.33 (rid). Ud refers to the degree of the polynomial form of the 
output polynomial, except step 1. That is, keeping the terms up to At”"* terms, 
where A is context-relevant coefficient. Ud = rA^. 

Definition 2.34 (n^p). refers to the degree of the polynomial form of the 
output polynomial at step 1. The reason for separation between and rid 
will be explained below. Uda is also the imposed polynomial degree for f{t). 

Definition 2.35 (polynomial coefficient error). Polynomial coefficient errors 
of the polynomial form of some g{t) are the errors caused by finite-precision 
limit placed when storing the digits of polynomial coefficients of each polyno¬ 
mial term Akt^ of the polynomial form of some g{t), with higher polynomial 
terms truncated. I will ignore this type of problem until the very end of this 
introduction section. 

Definition 2.36 {H{s), Hr{s)). For H{s), depending on context initially, it 
will refer to general linear filters, given in Laplace transform with variable s. 
However, for most of time iL(s) = l/(s -I- 1). iLr(s) = — l/(s — 1). 
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Definition 2.37 (step, step input). A step sp contains filter H{s), with the 
assigned initial condition 0 at t = Tsp.e (defined below), along with filter input 
and output. There are two steps used in this paper, and the second step takes 
the polynomial input from the polynomial output of the first step. The first 
step takes the polynomial form of f{t) as input. Step input does not refer to 
the heaviside step function input, but the input of some step sp. 

Definition 2.38 (T^p^e, n, r 2 , Vgp, a). Ti^e = ?'i, Tlj.e = ?'2 and Tgp^e = Tsp- 
a = . The rule for setting o, Vgp is given in Equation for all steps except 

sp = l,nd + 2,rid + 3. 

Definition 2.39 {parity{i)). parity{i) refers to even-odd parity of integer i. 

For odd integers, parity(i) = —1. For even integers, parity{i) = 1. Notice that 
the parity function is not 0/1-valued, but —1/1-valued. 

Definition 2.40 (high-frequency steady-state input(s)/response(s)). High-frequency 
steady-state input(s) and response(s) refers to non-zero angular frequency parts 
of some given steady-state input/response. (Recall that steady-state responses 
are defined in terms of Fourier series.) 

Definition 2.41 (plural form, part(s), input(s), output(s)). Recall that in the 
definition of high-frequency steady-state inputs, we used the plural form. We 
may have substituted the word “parts” equivalently. So why was the word 
“inputs” used? This is because by the linear time-invariant property of H{s), it 
is better to think in terms of a single-angular-frequency signal separately, when 
analyzing responses to steady-state input. This is why we use the plural term, 
so I wish this does not cause any confusion here. 


(r '|sp-i 

{sp-iy. 


n^-sp+l 

E 


2=0 


parityii) — 
z! 


1 


( 1 ) 


For step 1, a rule concerning a must be set. Thus, once a is set, this determines 
all Tgp, via Equation[^ except and Tn^+s- 


Definition 2.42 (step initial condition requirement). Step sp initial condi¬ 
tion (output-side) requirement is satisfied by the following: First, pass the 
polynomial form of step input to H(s), resulting in gsp(t), and get the value 
Wsp = gsp(Tsp). Expand to degree rid, and substitute Csp, and get value 
e ~ ^sp. Add the polynomial form of WspC f ^sp (to degree Ud) to gsp{t), 
which results in step output. 


Definition 2.43 (steady-state response, transient response). Steady-responses 
are the responses that can be written as Fourier series. Transient responses 
are the responses that cannot be written as Fourier series. These concepts 
are adapted from sinusoidal filtering literature, where transient response is the 
response that dies off to zero as t —^ oo. Here, transient response is defined 
more broadly - however, for purposes of this paper, these distinctions will not 
matter. A different definition is provided just for convenience of explanations 
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in this paper. In this paper, responses will refer to filter/step output, rather 
than input, unless explicitly noted, such as “input transient response.” Input 
transient response refers to transient response being used as filter/step input. 

Definition 2.44 (transient response/steady-state response restriction). In this 
paper, I will restrict transient response to those that may have chance of being 
generated by allowed graphs (transforming to f(t)). That is, there should be 
some steady-state response (including f(t) itself) input that generates transient 
response, and steady-state responses themselves must be allowed by the given 
input graph G. 

Definition 2.45 (steady-state input). This refers to steady-state response be¬ 
ing used as filter/step input. As said in the definition of steady-state response, 
if not noted otherwise with label “input,” steady-state response refers to fil¬ 
ter/step output. 

Definition 2.46 (transient response induced/caused/created (three words equiv¬ 
alently used) by steady-state input). This will refer to step output caused by 
steady-state input, as the name implies. If we talk about transient response at 
step sp induced by steady-state input at sp' where sp' < sp, this refers to the 
following. Step sp' steady-state input causes transient response Ae~'. At step 
sp' + I, suppose step input is Ae~' only. Then use the resulting step output 
of step sp' + 1 (which is transient response) for step input of step sp' + 2, and 
continue on until step sp. It is this step sp output (assuming Ae~' step input at 
step sp' -I- I) that we call transient response at step sp induced by steady-state 
input at step sp' . That is, transient response induced by steady-state input ne¬ 
glects the effects of steady-state responses all together from step sp' , and deals 
only with what happens by transient response of steady-state input at step sp', 
which pipe into step sp. 

Definition 2.47 (Integrative nature of HR filters). Consider H{s) = l/(s -I- 1) 
H{s) can be dissected into: H{s) = I/s—I/s^-|-l/s^ —l/s"*-!-... Notice how i7(s) 
is essentially a sum of cascades of integrators. This means that one does not 
have to worry about polynomial coefficient errors, as long as there is no input 
polynomial coefficient error, when one does not consider the initial condition 
change of the step (defined above). 

Now let us think about filtering f{t) so that one may obtain zero-frequency 
amplitude nu. Naively just using H{s) would result in huge error. This is 
because, angular frequency w is very large that finite-degree truncation error 
will be large. Relying on samples (and digital filters ) is not feasible, as they 
basically have the same problem - one can see this clearly by considering the 
polynomial interpolation of these samples. 

Thus, a clever mechanism is required, and this paper intends to provide that 
one. The integrative nature of HR filters would prove to be very useful. 
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2.2.10 CC2: clever mechanism 


Let us evaluate how a must be set, and why Equationwas set. Because of the 
flow of the discussion, in case any confusion arises regarding n^,!, please refer 
to sub-subsection [^.171 

Equation [2 is the truncated version of the following equation: 


a 




(sp- 1 )! 


= 1 


( 2 ) 


This comes from the idea that we would like to eliminate Ae * terms when the 
first step input is fco S K. Consider the first step output with the input fco: 


Oi(t) = ko — k^ae ‘ 


(3) 


where a is set by imposing the zero initial condition at t = ri. Pass this output 
to the second step: 

fco -I- 7 e“* — koate~* 

One would like to set 7 as zero. Separate oi{t) into fcg and —kQae~*. What 
Equation [ 2 ] says is that the Ae~* term of the second step output to —fcgae”* 
must be cancelled out by Ae~* term of the second step output to fcp. That is, 
at t = r 2 , one should make —koae~* = —k^ at t = r 2 so that it is cancelled 
out by fcoj which means the zero initial condition at t = r 2 is satisfied without 
creating additional transient responses. And so forth with other steps. 

So far, I have considered the case with the first step input only being ko = nh- 
But we also have high-frequency steady-state inputs/outputs, which need to be 
considered. 

The fundamental choice being made here is this: we will allow the first 
step transient response to high-frequency steady-state response to be 
extremely large coming mostly from “numerical errors.” That is, we 
will set ri and a to be relatively large. This by Equation implies a very 
small rsp for steps other than the first step, which also means there is very small 
numerical error. By numerical error, it refers to the error in calculating the value 
at t = Tsp of the step input passed to H{s) due to the finite (polynomial) degree 
restriction of the input. 

Now recall the Ud degree limit. Consider again input ko- This means that at 
step Ud + 1 , the output is: 

, koat'^‘^e~* 

ko -j- 

nd'- 

but here serves effectively as 1. Let us pass this output as an input to filter 
H{s) and we see that koodA'^e~^/nd\ disappears completely, and only ko remains 
the relevant input for the output. This is an important point, for we can define 
a cycle consisting of -I- 1 steps. 

After finishing the cycle, we have to re-define r„^+ 2 , as Equation no longer 
applies. Before doing so, let us consider what happens at the input of step 
nd + 2. Let us define r^, r^+i: 
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Definition 2.48 (r^, r^+i, fi). = r „^+2 = r(„^+ 2 ), r^+i = r„_j+ 3 . fi = nd+2. 
Assume that step fi input can be approximated by the following: 

AA i-iAp-* 

v{t) = kQ + ^ —-— = fco + Ha (4) 

i=0 

We will return back to the question of whether the approximation is useful for 
our analysis. 

Let us again divide the input approximation into two inputs: ko and Pia- First 
consider ko. The output by passing ko to step /r is: 

ko — kofie~*' 


which means /3 = e’''', or more correctly approximately with the right choice of 
and rid- Pass this to step /r + 1 and specify the desired output: 

fco — fco/3te“* 

Notice again that intention is eliminating term. This is done by: 

p^M + 1 

^-= /3 = e"'' (5) 

V+i 

(with the right choice of rid and and r^+i.) 

Let us specify the step /i + 1 output, given the approximation input ^{t): 

ko ~ koPte * + lib + lie + lie; + lie (6) 


where 116 , 110 , 11^,116 are defined by: 

lib = —ki [tr — 1] (7) 

He = fci \tr (e*) — I — t] e“‘ (8) 

Hd = fci [tr {e^r) _ 1 ] r^+ie~* (9) 

He =-fci [tr (e"'''+i) - 1 - r^+i] e"‘ (10) 


and tr(e’'=p) refers to truncating e* to degree nd and substituting t = rgp, and 
similarly with tr{e*), which truncates e* to degree but without substitution. 
Let us provide a system of equations of the constant term (cq) and At term (cit, 
A = Cl) of step /r + 1 polynomial output: 


ko H 
-ko 


ki 


o^M +1 _ 


- /ci 




r^,+l ^ r^+i 

Solving the system of equations results in: 


II = Co 

= Cl 


ko = 


Co 


/3 - 


Cl 


e’'^‘+ir^+i - + r^+i 


r^+i 


gv+i 


ev+i _ 


( 11 ) 


( 12 ) 


rih = ko 


Co 


Cl 


(13) 


To derive Equation 13 from Equation HA this of course requires the condition 


that /3 3> 1. Double-approximative equation 13 captures the dominant terms in 
denominators in Equation |12[ 
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2.2.11 CC2: role of Equation |12| and 

Definition 2.49 (Filter-H-initial value). Recall the definition for a step. A 
step first takes filter input, pass it to H{s) and impose the initial condition. 
Let us eliminate the initial condition and think as if we passed the step input 
to filter H{s). Then filter-iJ-initial value of step sp refers to the step output 
without imposed initial condition at t = rgp. 


Before going on, let me emphasize that Equation 12 will not be used for 
solving fcp. After all, there is danger that inaccurate ko = nt may be obtained. 
So what is Equation [T^ a ctually giving us, and what exactly is the approximate 
input v{t) in Equation ® 

Let us come back to step 2. Suppose step 1 input is kxe~*. Then passing this 
to step 2 results in: 

—kxr2e~* + kxte~* 


assuming numerical error in computing filter-iJ-initial value due to the polyno¬ 
mial approximation is small enough. Notice that r2 < 1. For step 3, the output 
becomes: 


kx ( + r2r3 


— kxr2te * + 


Consider how r 2 r 3 and /2 compares. Essentially, (ra)^ « r2, as both 
T 2 ,'I ’3 ^ 1. Thus, r2r3 is much less than /2. Continue these processes 
until step rid + l. It is now seen that for step p, v{t) is created assuming 
that for step 2 < sp < p, Ae“* term of step sp output was affected 
only by the / {sp — 1)! term of step sp output. By “approximative” 

relation: r2 « for step d < rid + I, we can see that the actual step p 

input’s magnitude is (1 -I- pd)W{t)\, where pd < 

Thus now we have obtained the assumption: 2”“^, which is to be used 

with pd to figure the effects of the parts truncated from actual step p input by 
using v{t). 

System of equations is now given as: 


/co + (1 + Pd)ki - 

-fcn^-(l+AidOfci 


'■m+i 

V+1 




V+i 


Co 

Cl 


(14) 


Note the added pd and pd' ■ Solve the equation and one sees that the dominant 
term equation, which is Equation |13| remains valid. 

Now coming to the question of why Equation 13 is useful. Recall that we 
have not so far considered steady-state response at step p + I and 
transient responses caused by high-frequency steady-state inputs of 
step 2 and beyond. What Equation [T3| allows us is to check whether they do 
matter. 


It is very easy to see that they do not. Again, recall Equation 13 


Co 


/3 


Cl 

/3 


Uh = ko 
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The polynomial term of degree rid of the steady-state response of high-frequency 
w at step fi+1 has the following (with denominator only capturing the dominant 
term): 

^ -ZW 

(jj'n-d+2 

Thus, it is clear that we have the right degree n-d at step ^ -I- 1 (this indeed 
is somewhat tautology via definition, but worth emphasizing in numerical con¬ 
cerns). For our relevant constant term and At term, this is more than sufficient. 
As the contribution to cq is much less than 1, our calculation of rih is not af¬ 
fected. Similarly with ci also. I will revisit this again, along with transient 
response analysis associated with these steady-state responses. 

Now contributions of transient responses not considered. The dominant con¬ 
tribution comes from step-2 transient responses induced by step 2 steady-state 
input (or step 1 steady-state response/output, equivalently). But e''''2”‘^n"/(w)^ 
can easily be made to be insignificant by choosing the right and c in f{t) = 
y{ct) so that these transient responses cause less than 1 contribution in magni¬ 
tude to Co or Cl. 

This completes our analysis (it may seem incomplete here, but the above sub¬ 
subsections will show that it is indeed complete) except the issues involving step 
fi and step yi + l transient responses induced by step and y, + 1 steady-state 
inputs, and requirements will be listed below. 

2.2.12 Numerical concerns, requirements 

• nd,i » rid, ndo » ri and Ud < ri 

— so that 2"‘^n"r„^+i <C 1 requirement may be satisfied: this require¬ 
ment deals with ir{t) and Equationj^so that analysis based on them 
can show that post-step-1 transient responses induced by post-step-1 
steady-state inputs and steady-state inputs themselves do not matter 
for our calculation of = fcg. 

— so that we may be able to shrink Vgp for 2 < sp < nd + 1 to satisfy 
uJTsp <C 1/n” in magnitude. 

• e’'''2"‘'n"/(w)^ 1/n^" so that post-step-1 steady-state inputs have do 

not create transient responses affecting our value of cq and ci significantly. 

• /3 = 1 so that Equation becomes valid, and so that above items 

are satisfied without a problem. 

• Finally, <C Ud so that filter-iJ-initial value error caused at step p is 
insignificant. 

2.2.13 Again, numerical concerns and other numerical concerns 

So far, we have not discussed how insignificant filter-i/-initial-value errors are, 
given our Vsp, ridp and rid- This is because they truly are very insignificant. 
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The below addresses that issue, plus provides a detailed summary of what was 
presented above. 

• Step 1 hlter-i/-initial-value errors do not matter. This is because ki al¬ 
ready incorporates these errors by assuming the actual A in the Ae~^ term 
at step 1. 

• Step 2 to step -I-1 to step + 1 

— Step 1 input fco always induces zero error: the rules for setting rgp 
ensure this. 

— We will not consider the fact that a cycle olud + l steps is formed due 
to polynomial truncation itself as being an error. What we consider 
as filter-i7-initial-value errors are mainly related to term trunca¬ 
tions. The exception is when we used e’’'' instead of the truncated 
one. 

— Let LOTsp ^ 1/n" for 2 < sp < fj,. Let and let r 2 = n?. 

Then, filter-iJ-initial value errors magnitude upper bound at step 
pL -\-l can be considered as -b 

for sum of post-step-1 transient responses induced by post-step-1- 
until-step-nd -I- 1 high-frequency steady-state inputs. 

— Steady-state responses have zero error, via the integrative nature of 
H{s). (HR filter) 

— The filter-iJ-initial value errors magnitude upper bound at step p 
and pL + 1 induced by step p high-frequency steady-state inputs are 
meaningless. Instead, we rather would have to consider in terms of 
actual calculated steady-state value at t = at step pL. 

— What was discussed in the above item applies for the filter-iL-initial 
value errors magnitude upper bound at step pi + 1 induced by step 
pi-\- 1 high-frequency steady-state inputs. 

— To summarize the condition above and below, e’’^‘2"‘^n"/u;^ ^ 1 for 
high-frequency w. 

2.2.14 High-frequency steady-state input/output at step pi and /i-fl 
and associated transient responses 

When LOTsp ^ 1 is not satished, as is the case with step pi, it is possible that filter- 
/i-initial-value is extremely large for high-frequency |w| 3> 1. Fortunately, this 
is not the case. Recall the following polynomial term of degree k (denominator 
only capturing the dominant term) for the step pi high-frequency steady-state 
output: 

, ,k±k ±k 

A^^^— =A—— 

^nd+2 

where |bl| < n". Simplifying analysis, one may form the upper magnitude 
bound of the calculated polynomial high-frequency steady-state output as (with 
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appropriate with /{r^)\ « n’''': 

ndrP>^nP 


We would like to make the above satisfy: 

ndrP>^e''f^rU 


< 1 




(15) 


(16) 


But this is satisfied by the constraint/requirement related to step-2 steady-state- 
input-induced transient response: 






< 1 


(17) 


As long as the above equation is satisfied, step ^ -I- 1 steady-state input also is 


insignificant for our calculation of Uh- Equation 16 can also be used for step 
/r-|-1, with transient response at step /r being used as input for step gL-\-\. This 
is because r^+i ^ 1. 

2.2.15 CC2: implementation values 


Let us refer back to the two lists (List 2.2.12 and 2.2.13) above. Let ri d = ^ 

. Let us satisfy Equation[l7|first (also the last item in List 2.2.13). Let us 
choose c = ^ in f{t) = y{ct). This gives us angular frequency range of '> 

to positive or negative angular frequency directions, plus of course 

zero frequency, which hosts our hamiltonian path angular frequency. With this 
chosen, let us choose ri, which by reference to rid also allow us to restrict rid^i- 
Let ri = which is set with consideration that l/(n^"' 1)^/" = \ jrY^ \ 

which gives us the range of rep for 2 < sp < y. This choice of ri satishes 
the requirement that 2”‘^n"r„j+i <C 1 (the fifth item, or the third sub-item 
connected to the second item, of List |2.2.13| and t he second item, or the Hrst 
sub-item connected to the first item, of List 2.2.12). Together, these satisfy the 


third and fourth item of List 12.2.121 The fifth item and the sixth item of List 


2 .2.12 was satisfied via the choice of rid and r^. 

This leaves us with the choice of Let us choose rid^i = 
Let us summarize: 


ridp 

II 

) O 

rid 

= 

ri 

= 


= 

fit) 

II 


(18) 


We are really done in analysis. 


2.2.16 CC2: extraction of rih 

Note that while Equation |12| was useful in showing significant contributions that 
must be used to calculate = rih, it nevertheless does not by itself accurately 
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compute Uh- 

But one can easily see that the way to compute rih then is as follows: have 
as (pseudo-)step two input and run all the steps using all the defined Vgp and 
imposed rid- 

Definition 2.50 (pseudo-step). Pseudo-step is defined as to distinguish with 
actual step operations. Pseudo-step does not have pseudo-step 1: pseudo-steps 
start with pseudo-step 2. The input is always at pseudo-step 2. Pseudo-step 
^1 + 1 output is used to determine (/ipp and (pi^i below. Thus, pseudo-steps are 
only used to determine coefficients in the final system of equations used to solve 
rih = ko- 

This gives us a complete picture on the relationship of kQ(pQ Q + = cq 

and ki(j)i^o + = ci. (^o,o and was already known (and equivalent 

as in Equation [m so the problems that get resolved here is on and 
Algorithm 1^ provides the implementation. In the implementation ki is replaced 
with zi to give the system of equations: 

ko4>o,o + = Co 

ko4>i,o + = Cl 

where Zi refers to step 2 transient response input created by high-frequency 
steady-state inputs at step 1. 

2.2.17 Revisiting Ud.i 

From the discussion, it may not be clear why must be different from 
This all comes from the fact that when is expanded and substituted with ri, 
adding higher-polynomial-degree polynomial terms of makes the calculated 
value diverge from actual value before reaching a certain degree. And it is 
the calculated that is used to determine A in Ae“* to be added to the filter 
output, given step input. When nd,i is around ri or less than ri, then calculated 
e~'^^ value is far off from actual 6“’'^ value. And because of this explosion, we see 
that big ri, given insufficient rid actually shrinks a in Equation This means 
r2 and other r^pS are affected, which now need to be much greater, causing large 
deviation in steady-state value (filter-i7-initial value error). 

This is why we set to be much greater than n^. 

2.2.18 For both Core Concepts: finite coefficient precision issue and 
implementation value of pi,p2 

So far, I have discussed the effects caused by the polynomial form of functions 
truncating polynomial terms of degree greater than rid- But polynomial coeffi¬ 
cients used for functions also need to be stored in finite and reasonable amount 
of digits. 

The idea here is incredibly simple. Recall our system of equations for cq and ci 
above. The equations say that as long as polynomial coefficient precision is rela¬ 
tively big enough, additional factors/deviations in cq and ci are too insignificant 
to values of fco and ki. 
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Definition 2.51 (polynomial coefficient precision: pi,p 2 ). Polynomial coeffi¬ 
cient precision digits. These are set to be more than sufficient. pi = binary 
digits (precision imposed on f{t) and grid intermediate calculations), p 2 = 
binary digits (Extraction/filter interemediate and final results). 

2.3 Algorithmic approach conclusion and time complexity 
analysis 

2.3.1 Core Concept Summary 

• A given undirected graph G = {V,E) was converted to f{t). Then, we 
formulated the signal processing problem as zero-frequency amplitude ex¬ 
traction of the Fourier series f(t). This zero-frequency amplitude gives 
one n^, the number of hamiltonian paths. 

• We directly work on the polynomial approximations of functions, instead 
of using samples. This makes us choose analog HR filters, as it provides 
precise polynomial coefficients for filter outputs, regardless of truncation 
degree chosen, as long as input polynomial coefficients are accurate. By 
the polynomial approximations, it refers to Taylor expansion around t = 0, 
truncated to degree This truncation degree is fixed for every filter 
output and input. 

• The main idea used in this paper is about eliminating Ae~* term when 
constant term was the input. That is, suppose f(t) = Uh- Then for every 
filter output, we would like to eliminate the effects of transient response 
at t = 0, which means eliminating Ae~^. But notice that when one passes 
rifi to filter H{s), transient response is inevitable. Thus we instead choose 
to eliminate the effects after one filtering operation (step) is completed, 
meaning Ae“* effects are eliminated at the second step. (Each step uses 
filter H(s) with atypical initial condition.) That is, we filter several times 
- which means implementing high-order filters as a cascade of first-order 
filters, but with atypical initial condition. (Typical initial condition is 
filter output being zero at t = 0. We instead choose to make filter output 
zero at t = r^p, where sp refers to step, or the ordinal of the filter in the 
cascade of filters used. 

• We of course have to consider high-frequency steady-state inputs of f{t). 
But eliminating high-frequency steady-state input-induced transient re¬ 
sponse effects at t = 0 is difficult (as we do not know exact w) and is not 
done. Instead, via clever use of Vgp and n^, we form a system of equations 
so that we may obtain nh- 

• What happens in Core Concept 2 is incredibly simple, even though details 
are not. Polynomial truncation at degree Ud means that output to input 
transient response e~*/nA disappears, meaning that -I- 1 steps form 
a cycle. Equationshows that if ri is very big, then r 2 is very small. Since 
we would like the polynomial truncation effects of steady-state responses 
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to be small for many steps as possible, this is the way we choose. This 
allows one to control transient response magnitude, except for the hrst 
step. 

• Step and /u + 1 work essentially in the same spirit as step 1 and 2, 
except that now we do not need to set /r to be very large. This is because 
steady-state response decayed significantly to allow us to take benefits 
of polynomial truncation. (We often think polynomial truncation only 
causes numerical errors, and these errors are in fact what this paper tried 
to overcome. But at step /x and fx + 1, they rather come as advantage in 
that we really do not wish high-frequency steady-state information and 
want them to be completely ignored as zero.) 

• Why is set is explained again in the separate sub-subsection and will 
not be mentioned again. The sub-subsection provides an insight on how 
Ud is set, along with other sub-subsections, and will be mentioned again. 

Another Core Concept 2 summary (a very short one without explanations): 

• The point basically is to set a system of equations that involves only 
n-h = ko and ki in the hrst-step-high-frequency steady-state input-induced 
transient response kie~*. This requires setting r^p, rid and nd,i so that 
other types of transient responses not caused by the first step steady-state 
input can be ignored. 

2.3.2 Time complexity analysis 

Even without going into deep analysis, it is clear by now that the above method 
provides a polynomial time algorithm for computing Uh- But for sake of com¬ 
pleteness, let us perform time complexity analysis. 

First, consider what time complexity would be like for filtering process. We used 
polynomial forms for hlter inputs and outputs. And in Laplace domain is 
1 /g/c-i-i^ which is a cascade of integrators. 

Thus, relative to the length of analysis we have done so far, actual hltering 
process is ridiculously simple. 

Let the temporary storage be the polynomial u(t), initialized with zero. Let p{t) 
be the polynomial form of inverse laplace transform of l/(s -I- 1) to degree Ud, 
and pi{t) be the polynomial form of inverse laplace transform of l/(s-|- 1 ) to de¬ 
gree ridp- For now let us consider normalized polynomial coefficients {B = Ak\ 
in At^). Some step polynomial term Bkt^/kl translates to simply moving the 
each normalized polynomial coefficient of the polynomial form p{t) by fc -I- 1 to 
the right and multiplying these normalized coefficients by B^. Then add the 
resulting coefficients to u{t). Continue the process from fc = 0 to /c = (or 
k = UdA in case of step 1) of step input. Then, get u{rsp), and add the poly¬ 
nomial form of —u{rsp)p(t)/p{rsp) to u(t). (pi{rsp) and pi{t) instead of p{rsp) 
and p{t) in case of step 1) 

Thus, for complexity purpose in O-time, time complexity is 0{nd(nd,i)^ip2)^), 
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with (p 2 )^ referring to multiplication of two p 2 -digit numbers, refers to 

the fact that there are at maximum polynomial terms that do integration, 
or polynomial coefficient shift, that have multiplicative Bk of the step input 
to be multiplied, and polynomial terms at maximum that the multi¬ 
plicative term is multiplied to. 0{nd + 3) refers to number of steps being used. 
Another nd,i refers to the number of multiplications for t^ terms (used when 
determining the value at t = Vsp), with k = nd,i at maximum (exponentiation). 
To calculate, 0(n^^°). Note that for simplicity, we actually derived time com¬ 
plexity that is much more than actual time complexity. In fact, the parameters 
used for this paper could have been less. 

Now consider the grid procedure. Time complexity is 0(n^ • (nd,i)^(pi)^), which 
refers to the dominant e®”* multiplication part of the inner loop of the grid. See 
the Algorithm for reference, n comes from the outer loop, the other n comes 
from the inner loop, excluding the inner loop of the inner loop (used for co¬ 
efficient additions), rid^i refers to the times/number of terms each polynomial 
coefficient of e™* needs to be multiplied to, the other 71^,1 refers to the number of 
such polynomial coefficients of e®"*. And (pi)^ of course refers to multiplication 
of two pi-digit numbers. This results in 0(n^°^) « 0(n^^°). 

So far, total time complexity was 0(n^^°), but note again that we have so far 
selected sufficient parameters, not efficient parameters. 

Now time complexity for determining r^p. Determination is done by Equation 
[2 which gives us a problem of solving the polynomial equation equalling to zero 
- that is find the root of the polynomial equation. There are at maximum 0(nd) 
polynomial terms, plus coefficients are 0(n^*^)-digits. Since there are many poly¬ 
nomial equation solvers that are of polynomial time to Ud and pi = (see, 
for example, will simply denote time complexity as PC{nd,Pi)- 

Now comes the extraction (of rih) part. But this part provide relatively insignif¬ 
icant time complexity. It was noted in sub-subsection 2.2. 16| that the equation 
and essentially the pseudocode we use is this: 


ko4'o,o + fci^o,i — Co 
ko4'i,o + = Cl 


( 20 ) 


(with of course, modification as in Equation[^ but this really does not matter.) 
And solving the system of equations is obviously pure multiplication and sub¬ 
traction and is largely insignificant. Rather, what is significant is determining 
(or more precisely, (pQ ^ in Equation 19) and (or more precisely, p'l j^). 
But this is simply running step 2 input e ^ until step p + 1, which gives us (pQ^i 
and Thus, total time complexity derived is 0(n^^°) -I- 0{ndPC{nd,P2)), 
which is polynomial-time complexity. 


2.3.3 Drawing the result from limited additional information 

It may seem at first that because there are exponentially many angular frequen¬ 
cies considered with equal angular frequency, there is no way one can efficiently 
draw out n^,. 

However, we do not need accurate information of amplitudes of f{t) at angular 
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frequencies other than zero. Furthermore, we know where maximum angular 
frequency with non-zero amplitude possibly lies, and where minimum non-zero 
angular frequency with non-zero amplitude lies. We also know that maximum 
sum of amplitudes and each amplitude can only be of magnitude n" at maxi¬ 
mum. 

Just like compressed sensing uses limited additional information - for example 
assumption of sparsity - to draw out results, limited information does allow us 
to get the desired result n^. 

2.4 Pseudocode for grid 

The pseudocode is provided by Algorithm [TJ 


Algorithm 1: Grid procedure for calculating f{t) 

Input: Graph G — (V, E), coeflicient precision pi — 

Output: The polynomial form of f{t) to degree rid within coefficient precision p 2 
n ^ \V\-, 

Let depth be d, V[l] be Zth vertex, wire be w; 

iu[Z][cZ] is a vector of polynomial coefficients for outgoing wires from depth d vertex V[l] to 

cZ + 1; 

d starts at 1, ends at n. At d — n 1, the algorithm stops; 

stands for the vector of polynomial coefficients up to polynomial degree nd\ 
e^^^’s polynomial coefficients are determined by Taylor expansion around t — 0; 

V[l] —— n\ where Z goes from 1 to n, with V matches to vertices of G; 
z * refers to multiplication of the polynomial z by the polynomial form of 
z * also truncates the multiplication result to polynomial degree rid', 

When calculating z * polynomial coefficient is calculated to pi digits; 

^ J2?=i yii'h 

initialize ii;[Z][l] — 
tZ ■«- 2; 

repeat 

repeat 

j ^ 1; 

repeat 

if iV[j],V[l]) G E then 

w[l][d] -(r- w[l][d] + w[j][d - 1]; 

3^3 + 1; 

until j > n; 

w[l][d] ^ w[l][d] * 

Z Z + 1; 

until Z > n; 

y(t) ^ (Er=i “[G"]) * 6“’“'**; 

fit) = 

(Z ■<— cZ + 1; 
until d > n; 


2.5 Pseudocode for multi-staged/multi-step filtering 

The algorithm is provided in Algorithmic The algorithm for Vsp is the polyno¬ 
mial root solving, presented in Equation ^ Normalized polynomial coefficient 
refers to A in At^/k\, and the definition of normalized polynomial form follows 
from the aforementioned. 
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Algorithm 2: Filter procedure 

Input: The polynomial form of /(f), coefficient precision p 2 — digits and Vsp 
Output: The polynomial form of o{t) with precision p 2 digits 
sp is simply s in this pseudocode with Vsp written as rs‘, 

All polynomial forms assumed to be of degree rid, except for step s = 1 which is of degree 

* refers to multiplication, with the product result truncated to precision p 2 digits; 

Let p{t) — e~*, and p[i] represents normalized polynomial coefficient to degree nd,i', 

Let temporary storage be j[] (vector), initialized with the normalized polynomial form of 

/(t); 

Let temporary storage be jp[], initialized with O’s; 

Let temporary storage be q and re, initialized with 0; 

Let m be initialized with nd,i', 

d, 2 ■<— 0; s ■<— 1; 

repeat 

if s 7 ^ 1 then 
m •«— rid’, 
repeat 

repeat 

if 2 + d+ l<m then 

ip[i + d + 1] •(- /p[i + d + 1] + j[d] * p[i]; 

until i > m; 

i <— 0; d <— d 1] 

until d > m; 

d ■«- 0: 

repeat 

w w jp[i] * {tsY/ i\] 
g -4— g + {PsY parity 
i <— i 1; 

until i > m; 

p i - w/q\ q^w,i -4— 0; 

repeat 

ipW <- /p[4 + M * parity(i); 
i ^ i+ 1; 
until i > m; 

2,d-4— 0; s-4— s + 1; 

/[] ^ ipD; ip[] <- 0; 
s s + 1; 
until s > rid + 3; 


2.6 Pseudocode for pseudo-steps 

The algorithm is in Algorithm]^ 

2.7 Pseudocode for extraction of Uh 

The algorithm is provided in Algorithm o{t) refers to o{t) in Algorithm 
(j)Q I and (j)i I used as in Algorithm 

3 Analog approach 

The analog approach is discussed completely separately from the algorithmic 
(digital) approach. Thus, both sections can be read in standalone ways. The 
analog approach is mentioned, because as seen in the algorithmic approach, 
analog signal processing was the foundation, instead of digital signal process- 
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Algorithm 3: pseudo-step procedure 

Input: None. Assumed input e“*. 

Output: 0Q and 

sp is simply s in this pseudocode with Tgp written as 
All polynomial forms assumed to be of degree rid', 

* refers to multiplication, with the product result truncated to precision p 2 digits; 

Let p(t) — and p[i] represents normalized polynomial coefficient to degree rid', 

Let temporary storage be j[] (vector), initialized with the normalized polynomial form of 

Let temporary storage be jp[], initialized with O’s; 

Let temporary storage be q, w, initialized with 0; 

Let m be initialized with n^; 
d, i — 0; s — 2; 

repeat 

repeat 

repeat 

if i + d+ l<m then 

jp[i + d + 1] -i- jpli + d + 1] + jld] * p[i]; 
i -4— 2 + 1; 
until i > m; 

1 -4— 0; d -4— d + 1; 
until d > m; 

d -4- 0; 
repeat 

re -4— ic + jp[i] * i'f'sY 
g -4— g* + {PsYparity 
■i -4— i + 1; 

until i > m; 

< - w/q,w,i i— 0; 

repeat 

ip W jp[i] + parity(i); 

2 -4— 2 + 1; 

until i > m; 

2,d-4— 0; s-4— s+1; 

i[] <- ipD; ipQ <- o; 

s -4— s + 1; 
until s > rid + 3; 

<io,i i[0]; e- i[l]; 


ing - basically, polynomial (analog signal processing foundation) vs. samples 
(digital signal processing) for the form of inputs and outputs in the algorithmic 
approach. Thus the analog approach is more “natural” way of understanding, 
even though it has less relevance in reality and may not even be operative. Any¬ 
one solely interested in the digital computer algorithm for obtaining the number 
of hamiltonian paths in undirected graph G = {V, E) will gain nothing from this 
section. Nevertheless because of the history behind the algorithmic approach, 
the analog approach is mentioned also. 

Before going on, in practice and in theory there is no good frequency multiplier 
except for a signal of single frequency sinusoid. Thus, the use of frequency mul¬ 
tiplier is a conceptual one. An alternative circuit formulation, used somewhat 
identically for the algorithmic approach, is also shown and used. 

In general, analog methods do suffer from several problems, and this should be 
in consideration when reading off the results in this section. 

Graph G = {V, E) is defined with a set V of vertices along with a set E of 
undirected edges connecting vertices. We will denote a walk by the following 
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Algorithm 4: rih extraction procedure 

Input: j, o(t), \V\ = n, r^, 

Output: rifi 

Use system-of-two-equations-solver for Equation |19[ with o[0] = cq, o[1] = c\ and given 
inputs; 

Round the resulting to the nearest integer; 


formalism: a — b — c where a, 6, c are vertices and — represents edges. Generally, 
a, &, c will be represented with positive integers. 

n = riy is the cardinality of V, rig is the cardinality of E. A n-walk is dehned to 
be a walk with n vertices. This is the only class of walks we will have interests 
in this paper. 

We will re-interpret Hamiltonian path existence problem using a nx n grid and 
others. 

Definition 3.1. The grid contains n vertical columns. Each column contains 
n vertices, and the vertices in the same column are not connected by wires. 

Definition 3.2. All vertices are numbered with positive integers greater than 

1 . 

Definition 3.3. Each wire transmits a voltage signal f{t). For our consid¬ 
eration, location does not matter, so all of our signals are solely function of 
time. These signals can be transformed into the Fourier transform frequency 
representation. 

Definition 3.4. As part of lumped circuit assumption, we will assume that 
wires have no time delay, (ideal wire) 

Definition 3.5. For each vertex x at column a > 1, if vertex y satisfies {x, y) G 
E or {y, x) G E, vertex y frequency multiplier (or oscillator, in case of a — 1 = 1) 
at column a — 1 is connected by a wire to the sum operator at vertex x /column 
a. 

Definition 3.6. As we allow self-loops, while {x,x) ^ E, vertex x at column 
a — 1 is connected by a wire to the sum operator at vertex x/column a. 

Definition 3.7. Each vertex x at column 1, the first column, only has an ideal 
oscillator, transmitting to wires connected to the second column. 

Definition 3.8. A sum operator just sums up the signals transmitted by wires. 

Definition 3.9. Each sum operator at vertex x/column a is connected to a 
frequency multiplier at the same vertex/column, with frequency multiplication 
factor of X. Frequency multiplier transforms -I- -I- ... into -|- 

^ixW2t _|_ 

Definition 3.10. At column n, after signals pass through frequency multipliers 
connected to sum operators, any wire incident from column n is connected to a 
final sum operator, which produces the final signal y(t). 





Signal Processing Understanding of Hamiltonian Paths 


27 


Thus it is clear that we need n{n — 1) + 1 sum operators (or adders, equiv¬ 
alently) and n(n — 1) frequency mnltipliers for the grid above. The number 
of wires are dependent on E, but the maximum number of wires required is 
n^(n — 1) -b n{n — 1) -b n, where the last n comes the wires that connect column 
n to the last sum operator, and n(n — 1) comes from the wires that connect a 
single sum operator to a single frequency multipliers. 

The output of the circuit grid defined above is y{t), as mentioned above. Let 

V = {vi,V2, ■.,Vn}. 

Definition 3.11. The final sum operator, which produces the signal i/(t) is 
connected to the ideal mixer M, which outputs the product of i/(t) with 
where u = viV2V3..Vn- In Fourier transform, this is equivalent to converting 
V(uj) with F(a; -b u), where V(uj) is Fourier transform of y(t). Let the output 
of M be k{t). 

From the above, it is clear that Ce™* inside y{t) represents hamiltonian 
paths, with C representing the number of hamiltonian paths. In k{t), frequency 
0 represents hamiltonian paths, as all frequencies are shifted left by u. 

Because our chosen low-pass filter will be first-order, we will also pass k{t) to a 
frequency mnltiplier that multiplies frequencies by Vn^'^ where is the greatest- 
numbered vertex, to ensure that the frequencies other than zero frequency parts 
of k{t) will be sufficiently high frequencies. (Multiplying zero by Vn'^'^ is zero) For 
higher-order filters, like third-order filter, this additional frequency-multiplying 
process will not be needed. We will call the resulting signal j{t). 

As a side note, instead of having input tape in Turing machine, we have to 
re-wire n x n grid every time graph input changes. This n x n grid serves as an 
input to the system involving a low-pass filter. 

3.1 Restriction on vertex indices 

However, a close look will reveal that it is required for us to restrict on vertex 
indices. Hamiltonian u may be decomposed into a product of n numbers that are 
in V, and yet all these numbers may not be distinct, required for u to represent 
hamiltonian paths. One simple way to address this problem is by required all 
vertex indices to be prime numbers. For simplification, assume that vi = 2 and 
Vn = Pn where pk represents fcth prime number with pi = 2. It is known that 
Pn < n(lnn -b Inlnn), shown in Rosser (1941). Thus we only need to check 
non-exponential number of natural numbers to obtain n prime numbers to be 
used as indices for vertices. 

3.2 Analog approach: low-pass filter 

Now that we defined the final output k{t), the question is how we process k{t) 
to give us some information about the number of hamiltonian paths, or C. 
To do this, we pass it to a low-pass filter. But we cannot simply assume an 
ideal low-pass filter, represented by H{uj) = rect{oj), where rect{uj) = 1 for 
—0.5 < uj < 0.5 and rect{u}) = 0 otherwise, because there is no such an ideal 
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Vin 



Miut 


Figure 3: A first-order RC low-pass filter 


filter even to the approximate level. 

Thus we will choose a simple physical first-order RC low-pass filter, described 
in figure]^ (The figure is in public domain.) 

By Kirchhoff’s Voltage Law, the low-pass filter in figure 1 has the ODE of: 


dUout 

dt 




T 


where r = RC. As this ODE is linear, to figure out the behavior of this low-pass 
filter, we first consider Vn = De**^*, where oj is some arbitrary frequency. 

Using initial capacitor voltage condition at the starting time t = 0 as Vout-.t=o = 

0 , 


^out 


D 

1 -I- iWT 



Assume that > n -|- 1. Also for calculation convenience, assume that r = 
RC = 1. In steady state t = oo, because every w of j{t) except zero is greater 
than/equal to and the total number of walks in C with n total vertices 

can only have maximum of n” n-walks, j(oo)’s value mostly comes from the 
hamiltonian/zero-frequency part. Other frequency parts only contribute less 
than 1/n^" in magnitude. Thus at time oo, the number of hamiltonian paths is 
discovered from the magnitude of j(oo), |j(oo)|. However, calculations must be 
done on finite time, so the steady-state case only forms a background for our 
discussions, not the main part. 

Note that in ordinary signal processing, keeping phase errors small is very im¬ 
portant, but for the use of signal processing tools to analyze hamiltonian paths, 
phase errors are not of any concern. 


3.3 Time Complexity of the Circuit 

But moving to the finite time is simple: figure out the time when decays 

to . Then high frequency parts only contribute a negligible value to j{t). 

Now since r = 1 assumption is made, set equality 6“*“= = 1/n"^". Taking the 
natural log to each side, tc = 4nlnn < 4n^. Thus, the critical time, which is 
when the exponential decaying factor decays to l/n'^", increases approximately 
linearly as the size of input n increases. 

After this critical value, the value of \j{t)\ can simply be sampled by a digital 
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computer to get the number of hamiltonian paths. Note that theoretically 
only one sample is required to measure the number of hamiltonian paths. This 
is because frequency 0 does not have any oscillating part, and thus will have 
constant offset relative to 0. 

Thus time complexity of the circuit to solve the number of hamiltonian paths 
is 0(n log n), which is smaller than 0{n^). 

3.4 Size and Time Complexity 

The above demonstrates that the number of needed components and needed 
time does not grow exponentially as the input graph size increase. All the values 
used in the circuit does not require exponentially-growing number of digits in a 
digital computer, as the graph size increases. 


3.5 Alternative circuit formation 


In this section, I will describe another way of building a circuit that represents 
a graph. This method eliminates the use of frequency multipliers, and replaces 
them with ordinary multipliers. 

Start with the original idea that each vertex x at column 1 transmits to the 
wires x at column 1 are connected to. All wires going to vertex y at column i > 1 
are first met with a sum operator, but now followed by an ordinary multiplier 
of sum X The method will be explained in detail below. 

Definition 3.1 3.2[ 3.3 3.4 will be used as before. Section [O] no longer applies 
and is replaced with the restrictions imposed by the following definitions in this 
subsection: 


Definition 3.12 (The set V of vertex numbers). The set V is defined as F = 
{n, ..., n"}, which represents the set of vertex numbers (or equivalently vertex 

indices), with \V\ = n, the number of vertices. 

Definition 3.13 (n-walk). A n-walk^ = (^i,'C 2 , ^n) with G V and (^i, ^i+i) G 

E or = ^i+i, a list, is a walk that has n vertices. A n-walk may contain self¬ 
loops or loops. One may consider a n-walk as a list of n vertex numbers that 
may contain one vertex number more than once. 

Definition 3.14 (Permutation of a list). A permutation of a list is a re-ordering 
of list elements of 

Definition 3.15 (Uniqueness of n-walk frequency). Let a n-walk ^ be ^ = 
(?i!C 2 i which is a list. Let ui = is a unique n-walk frequency 

of G if it can only be the sum of some permutations of one list. 

Lemma 3.1. For V = {n, n^,..., n"}, there cannot exist a n-walk frequency 
such that it is not unique. 

Proof. The proof is simply the basis representation theorem, except that the 
case where n vertex numbers that are same are in the list. In such a case. 
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u! = n ■ n\ But then uj = = 1 • and ^ = (n®+^) is the only possible 

alternative representation of uj. But the alternative list only has one vertex. 
Thus, there cannot exist a n-walk frequency that is not unique. □ 


Definition 3.5 needs to change as follows: 


Definition 3.16. For each vertex x at column a > 2, if vertex y satisfies 
{x,y) £ E or {y,x) € E, vertex y mixer at column a — 1, which multiplies 
to a signal it receives, is connected by a wire to the sum operator at vertex 
x/column a. In case of each vertex x at column a = 2, if vertex y satisfies 
{x,y) £ E or {y,x) £ E, vertex y oscillator (output of e®^‘) at column 1 is 
connected by a wire to the sum operator at vertex x/column 2. 

Definition |3.6[ |3.7| and |3.8| are kept. Definition |3.9| and |3.10| change to the 
following: 


Definition 3.17. Each sum operator at vertex a;/column a > 2 is connected 
to a mixer at the same column and the same vertex, which shifts frequency by 
X. A mixer, with shift factor of x, transforms e®®"i‘ + e®®"^* + ... into e®(®®'i+^h _|_ 
e®(®" 2 +®)* + ..., because it multiplies e®“* to the signal it receives. 


Definition 3.18. At column n, after signals pass through mixers connected to 
sum operators, any wire incident from column n is connected to a final sum 
operator instead, which produces the final signal y{t). 

Complexity remains the same: one needs n{n — 1) + 1 sum operators and 
n{n — 1) mixers/multipliers, (multipliers here are not frequency multipliers, but 
ordinary signal multipliers) The number of wires required remains the same. 
Definition |3.11| changes to the following: 

Definition 3.19. The final sum operator, which produces the signal y{t) is 
connected to the ideal mixer M, which outputs the product of y{t) with e”®®®* 
where u = Vi + V 2 + + .. + Vn, with Vi £ V. In Fourier transform, this is 

equivalent to converting Y(uj) with Y{uj + u), where Y{uj) is Fourier transform 
of y{t). Let the output of M be k{t). 

Now k(t) has zero frequency as its hamiltonian path frequency, as in the 
original formulation. 

One may choose to add frequency multiplier after the final mixer M so that 
a simple first-order low-pass filter can be used. However, one may instead 
choose to increase the difference between each vertex number, such as H = 
{n, n®®, n^®®,..., n®® }. This way, one does not have to add an extra frequency 
multiplier, which is likely to diverge from its ideal behavior, as I will discuss. 


3.6 Real deviations: perils of high frequency and frequency 
multipliers 

While the system described above is a physical system, not just a logical sys¬ 
tem, it is nevertheless still an ideal physical system. Oscillators are not perfect 
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oscillators, resistors and capacitors are not ideal ones, wires have impedance. 
Thermal effects may change system properties. 

But the most fundamental problem is the fact that the systems above are based 
on lumped-circuit analysis. Lumped circuit analysis works for low frequencies, 
because the length of wires can be made short enough to satisfy lumped-circuit 
assumptions. But one cannot shorten wires forever, and this makes lumped- 
circuit analysis to break for high frequencies. No longer discussion of lumped 
capacitors, resistors and inductors becomes a simple one. 

Many problems, whether small or not, require more details and are left out here. 
Future papers will address these issues. 

4 Conclusions 

For the algorithmic/digital approach, the conclusions are provided in the sepa¬ 
rate subsection. Note, or recall, again that digital and analog approaches were 
discussed separately and thus they can be read separately without affecting 
understanding. Thus, if the interest of the reader is solely on the algorithmic 
approach, one may skip reading the entire analog approach section and this 
conclusion. For the analog approach, the conclusion is that while ideal analog 
computing models - one like GPAC - provide a useful tool, they nevertheless 
fail for the problem we would like to address for all realistic implementations of 
analog computers. The question thus is how close can a real analog computer 
be to ideal computing models (this involves sciences, like physics, chemistry 
and biology) and under current realistic limitations, how improvements can be 
made. Some possible improvements have been provided in this paper, but rest 
are still left to future papers. Furthermore, while the digital approach is inspired 
by analog signal processing methods, using the digital approach may in the end 
be much better than the analog approach. 

4.1 An alternative model for analog computations 

As we have seen in the paper, analog signal processing is the foundation and 
is easily converted to an analog computer circuit. However, because of error 
build-up in analog computers, digital computing is usually preferred for large n 
(and in practice, no one uses an analog computer). 

So far, the GPAC ideal analog computer model was used for the analog ap¬ 
proach. However, one may also think of distributed analog computation models 
[5]. The idea is to provide a compromise between analog and digital approaches. 
The GPAC model assumes infinite precision digit, which is not the case in prac¬ 
tice. We may try to do get the best analog precision for wire precisions and 
components, but one may instead opt to distribute some of the digits. For 
example, if one requires 50-digit precision, then one may divide this precision 
into 10 5-digit-precision wires/components, with components interacting with 
each other (and processors are analog, instead of digital). As [5] shows, this 
provides improvements in analog signal-to-noise ratio. This, however, does not 
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completely eliminate significant noise accumulation over “stages” (to refer to the 
cases in this paper, “depths” or “levels” in a grid) as number of stages passed 
by increases. This makes us think of the hybrid approach in [5], which involves 
digital restoration and use of A/D/A. 

Here, one sees some common grounds with what was done in the algorithmic 
approach, even though analogy will not be perfect. I will leave more detailed 
arguments to future papers, and conclude this paper. 
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